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Abstract. We study the hadron-quark phase transition in a molecular dy- 
namics (MD) of quark degrees of freedom. The hadron state at low density and 
temperature, and the deconfmed quark state at high density and temperature 
are observed in our model. We investigate the equations of state and draw the 
phase-diagram at wide baryon density and temperature range. We also discuss 
the transport property, e.g. viscosity, of qq matter. It is found that the ratio 
of the shear viscosity to the entropy density is less than one for quark matter. 
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1. Introduction 



Recently experimental results which suggest properties of quark matter at high 
temperature have been obtained [E El E]- They tell us that the shear viscosity of 
quark gluon plasma (QGP) is very small and the phase is liquid-like. A theoretical 
study by quenched lattice QCD was performed on the matter viscosity at above 
the critical temperature [ 4 . The viscosity per entropy density 77/s was found to 
be much smaller than 1, which indicates the matter is extremely smooth. However, 
dynamics of quarks and anti-quarks, which seems very important, was not consid- 
ered in the study. In this paper we study the properties of QGP by employing a 
molecular dynamics (MD) which is a natural framework for describing both statical 
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and dynamical aspects of the phase transition We first show thermodynamical 
behavior of quark matter at finite densities. Then our preliminary results on the 
properties of matter at finite temperature and zero baryon density are presented. 

2. Framework 

We extend our previous MD model for zero temperature to finite temperature. 
We prepare a cubic box having a length of 4 - 10 fm, in which the same number 
of u and d quarks are contained, and impose a periodic boundary condition. To 
simulate quark many-body system at finite temperature, we give appropriate initial 
values of position and momentum of i-th quark, and P.;, and solve the following 
equations, 

where H is a Hamiltonian used in Ref. The equations of motion are solved for 
several thousands fm/c after the system reached equilibrium. 

We simulate two kinds of quark matter, i.e., a matter with finite baryon density 
and qq matter with zero baryon density. The first one does not include q even at 
finite temperature as a first trial. In the calculation of the second case, qq matter, 
we slightly modify the Hamiltonian, i.e. meson-exchange terms [|5] are switched off. 
The number of quarks and anti-quarks N q for qq matter is determined as follows: 
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where \i is a chemical potential, f2 a thermodynamical potential, T a temperature, 
V a volume and m denotes a quark mass of 300 MeV assumed in this calculation. 
Note that this formula is for a free Fermi gas and no effect of interaction is included. 
Therefore the number of quarks and anti-quarks is not consistent with the present 
MD simulation. 



3. Results and Discussion 

At first we discuss properties of quark matter. Figure 1 shows temperature de- 
pendence of the energy per baryon at 0.5pa and IpQ. The effective temperature is 
calculated from the momentum distribution [|S] as, 

1 JL,PidH 

i— 1 

In our MD simulation, quarks form baryons at low temperatures. At a certain 
temperature a baryon-quark transition occurs and quarks are released from baryons. 
This can be seen as the change of specific heat in Fig. 1. We define the critical 
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Fig. 1. Temperature dependence 
of energy per baryon for 0.5po (thin 
line) and lpo (dark line). The 
arrows indicate maximum-specific- 
heat points. 



Fig. 2. Phase diagram on the 
density-temperature plane. The 
saturation density of the symmet- 
ric nuclear matter is indicated by 
Po- 



temperatures by maximum-specific-heat points indicated by arrows. The critical 
temperature becomes smaller as the density increases. We draw a phase diagram 
of quark matter in Fig. 2. Baryon phase is realized at the lower side of density 
and temperature, and quark phase at the higher side. Note that qq pair creation 
and annihilation are neglected in this model. If thermally excited qq pairs were 
considered in the model, the boundary of baryon and quark phases might move 
to the lower side of density and temperature: when the space is filled with more 
quarks and anti-quarks, deconfinement of a hadron can occur more easily due to 
interaction with other quarks/anti-quarks in the environment. 

Let us next discuss thermodynamical properties of qq matter. The shear vis- 
cosity r\ can be calculated by using Green-Kubo relation 

V= T^vIo Jx v^' J ^ dt > (4) 
where J xy is the stress tensor in the x-y direction [Ej, 

To discuss the ratio of the shear viscosity to the entropy density, thermodynamic 
quantities in MD are calculated as 

e + P 1 / dH \ 

s = ^r< ™ =NT « + iiZ**\-m;)> (6) 

i—l x / 

where P is a pressure, s is a entropy density and e is the energy density. Simulation 
is performed by using a box size of 10 fm for low temperature region and 4 fm for 
high temperature region. Figure 3 shows the ratio of shear viscosity to the entropy 
density. In the low temperature region, where quarks form mesons, the ratio rj/s 
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Fig. 3. Temperature dependence 
of the ratio of the shear viscosity to 
the entropy density for qq matter. 
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Fig. 4. Radial distribution func- 
tion of qq matter at 120, 140, 160 
and 180 MeV from the top. 



is larger than one. The value decreases as the temperature increases. Around 130 
MeV, where quarks are deconfmed, rj/s takes a value smaller than one. This value 
is of the same order as other studies [0J IB]- Above 130 MeV, rj/s increases linearly. 
This means that the motion of quarks are more dominant than the interaction at 
high temperature. If an effect of the asymptotic freedom were taken into account, 
the ratio rj/s at high temperature would get larger. 

Radial distribution functions of relative distance of two quarks at 120, 140, 
160 and 180 MeV are shown in Fig. 4. They are normalized to one in the case of 
completely uniform distribution. Bumps near zero distance are due to the meson- 
like correlation. Change from the meson phase to the quark phase is seen as the 
temperature increases. However, since the number of quarks and anti-quarks N q is 
not consistently introduced in the present model, the critical temperature or other 
quantities of phase transition cannot be discussed quantitatively. 

The results of this paper are still on a preliminary level and longer simulation 
time and larger box size are necessary for more quantitative discussions. Investiga- 
tion of viscosity for matter of finite baryon density is in progress. 

Y. A. is grateful to H. Kaburaki and T. Hatsuda for their valuable discussions. 
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